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Abstract 

Measurements of energy separations in highly charged ions can in many 
cases nowadays be performed with very high accuracy — an accuracy that 
sometimes cannot be matched by the corresponding theoretical calcula- 
tions. Furthermore, it has recently been demonstrated (Chantlor et al. 
Phys. Rev. Lett. 109, 153001, 2012) that there is a systematic deviation 
between experimental and theoretical results for the X-ray lines of 
medium-heavy heliumlike ions. We have during a number of years been 
developing a general procedure for energy-dependent perturbative calcu- 
lations, which opens up a unique possibility of incorporating the energy- 
dependent QED perturbations into the all-order many-body perturbation 
expansion in a rigorous way. Such an expansion will yield several impor- 
tant effects, never before accounted for in this type of analysis, which is 
expected to increase the theoretical accuracy considerably. Calculation of 
some of these effects have been performed at our laboratory in Gothen- 
burg, and numerical results are given. Further work along this line is now 
in progress. To what extent the improved procedure might explain the 
discrepancy found by Chantler et al. remains to be seen. 

1 Introduction 

The theory of quantum electrodynamics (QED) has been very well tested for 
simple systems, like the free electron (anomalous magnetic moment) or the 
hydrogen atom (Lamb shift). For multi-electron systems, on the other hand, 
the corresponding tests are much less complete. There is presently a great 
research interest in testing QED, particularly of highly charged ions, which can 
be produced in good intensities at large accelerators, Uke that at the GSI facility 
in Darmstadt. 

In recent years very accurate experimental data have been produced for 
a number of highly charged few-electron ions, as well as for some inner-shell 
transitions of many-electron systems. The corresponding theoretical analysis. 
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however, can in some cases not match the experimental accuracy. In Tables ([T]) 
and ([2]) we show two examples of this situation, a transition in He-like silicon, 
measured with laser techniques by Myers et al. [T] , and the K„ X-ray lines in the 
copper atom, studied by Deslattes et al. j2] The theoretical estimates are made 
by means of two-photon exchange (Artemyev ) or relativistic many-body tech- 
nique with analytical QED energies added (Plante j4j, Chantler/Grant |5J). The 
experimental results are here 1-2 orders of magnitude more accurate. 

Furthermore, Chantler et al. [E* have recently found ''^Evidence for a Z- 
dependent Divergence Between Experiment and Calculation'^ for the lines in 
medium-heavy ions of several times the combined experimental and theoretical 
uncertainties (see Fig. [T]). The result for He- like Ti, studied by Chantler et 
al. at the EBIT facility at NIST is shown in Table [3] Here, the difference 
between theory and experiment is three times the experimental uncertainty. The 
theoretical results are obtained by Artemyev et al. [Sj and by Plante et al. [3], 
and the former claim much higher accuracy than the experimental result. 

The calculations of Artemyev et al. include two-photon exchange in a com- 
prehensive way but do not contain electron correlation beyond second order. 
The calculations of Plante et al. contain all-order electron correlation effects 
but are less complete in the QED part. 

It is evident that it would be possible to test QED effects on a deeper level, 
were more accurate theoretical results available. We have recently developed 
a procedure for energy-dependent perturbation theory [II [3 El [TOl HU [121 US] ; 
which for the first time will make it possible to include QED perturbations in an 
all-order perturbation expansion in a rigorous way. This will take account of new 
effects beyond two-photon exchange that have previously not been considered. 
This is expected to improve the theoretical accuracy considerably. 

Calculations of some higher-order QED-MBPT effects of the type mentioned 
above have been performed at our laboratory, so far restricted to the ground 
state of He-like ions. In his Ph.D. thesis work Daniel Hedendahl [S] calculated for 
the first time so-called non-radiative QED interactions (retardation and virtual 
pairs, see Fig. |4|, combined with all-order electron correlation (see Table |4] 
below). The inclusion of radiative QED effects (self-energy, vertex correction, 
and vacuum polarization) into the all-order perturbative expansion is presently 
under way at our laboratory. 

The outline of the paper is as follows. In the next section we shall sum- 
marize the well-known many-body perturbation theory for energy-independent 
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Table 2: Ka X-ray data for copper (in eV) 
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Table 3: The transition Is^ ^S^ - ls2piPi for He-like Ti (in eV) 
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perturbations, which forms the basis for the project. In the following section 
we shall briefly describe the newly developed theory for energy- dependent per- 
turbations, which is the instrument for incorporating the QED effects into the 
all-order perturbation expansion in a rigorous manner. In the final section we 
shall indicate how this can be accomplished and give some numerical results. 



2 Standard MBPT 
2.1 Non-relativistic MBPT 

As a background we begin by summarizing the standard time-independent per- 
turbation theory (see, for instance 14J). We consider a number of solutions to 
the non-relativistic Schrodinger equation, known as "target states", 

H-^'^ = E"^"' , (1) 

where the Hamiltonian is given by 

N 

H = Y,hs{i) + Vc. (2) 

1=1 

Here, hg \s the Schrodinger single-electron Hamiltonian and Vq is the electo- 
static interaction between the electrons 

For each target state there exists a zeroth-order or model function, 'I'g . A 
wave operator, VL, transforms the latter to the exact solutions 

= (4) 
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Figure 1: (a) The experimental data for He-like Ti. (b) The difference between experimental 
and theoretical Kq X lines for medium- heavy helium like ions). (From Chantler et al. [6], in 
eV.) 



The model functions form a model space, and in intermediate normalization 
these functions are the projections of the exact solutions on this space 



(5) 



We partition the Hamiltonian in the standard way into a model Hamiltonian 
and a perturbation 

H = Hq + V. (6) 



The wave operator satisfies the generalized Block equation 



n,Ho]P = Q{vn-fiw)p, 



where W is the effective interaction 

w = Ves = pvnp. 



(7) 



(8) 



This form of the Bloch equation is valid also when the model space contains 
different energy levels. 

In the general case we can separate the Bloch equation into one equation for 
each energy level £ of the model space 



np£^rQ{£){vri-nw)P£, 



(9) 
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Figure 2: Two-electron ladder diagram (left). In the corresponding counterterm (right) the 
energy parameter of the wave-operator part is changed to the energy of the intermediate 
model-space state {€'). The vertical lines represent electron orbitals/propagators and the 
dached lines the electrostatic interaction. Time is supposed to flow upwards. 



where Tq{8) — Qr{£), Q is the projection operator for the space outside the 
model space and 



10(^1 



(10) 



is the resolvent with summation over aU eigenstates \i) of Hq with eigenvalue 
. Expanding the equation ^ order by order, yields 



■m) 



(11) 



The last term in the Bloch equation (|9| and in the expansion ( 11 ) is usually 
referred to as " folded" , since it is conventionally represented by folded diagrams. 
It can also be regarded as the finite remainder, when the singularities — due to 
intermediate model-space states — of a ladder expansion are eliminated, which 
we shall illustrate as follows. 

We assume we have a wave-operator in the form of a "ladder" expansion 
(see Fig. [2j left) 

T{£)VT{£)VT£)V (12) 

where one of the intermediate states lies in the model space with an energy £' 
{Pg'), almost degenerate with the state we start from with energy £. The ladder 
then becomes 

rQ{£)v---^^vrQ{£)v---Ps, (13) 



which can be expressed 



n^'^\£)-^^W^''-"''>{£)P£, 



(14) 



assuming there are n interactions in total and m after the degeneracy. Since we 
have assumed that the interactions are energy independent, the wave operator 
and effective interaction depend on the energy only through the resolvents. 
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The expressions above are quasi-singular, and to eliminate the singularity 
we add a counterterm, (see Fig. [2] right) 



(15) 

where the energy of the wave operator is modified to E' . But the first term in 
the Bloch expansion (fl| contains counterterms from previous orders, namely 



The difference then becomes 



or 



{Tq{E) - VQ{E'))vn'^'''-^\E')j^ VF("-")(£)P£ 
- TQ{E)n^'"'\E) W^("-")(£: )P£, 



(16) 

(17) 
(18) 



which is the folded term in the expansion (11 1. 

If we assume that the interactions can be different in each order of the 
expansion, Vi, V2, Vis, then it follows that fl^™) should be formed by the m 
last interactions and VF^""™) by the remaining ones. 



2.2 All-order expansion 

By separating the wave operator into one-, two-,... body parts, 

17 = 1 + f7i + 1^2 + • • • , (19) 

the Bloch equation ([7| can be separated into a number of coupled equations, 

[^l,,,Ha]P = Q{V^l-nw)^P. (20) 

Normally, the two-body part dominates heavily, and this part leads to the pair 
equation. Efficient numerical methods have been developed for solving this equa- 
tion for atomic systems |15l 116] . Solving this equation iteratively (neglecting 
the minor one-body part), leads to the all-order pair Junction^ illustrated in Fig. 
[Sj For a two-electron system, starting from hydrogenic orbitals with no electron 
core, this will ultimately lead to the exact non-relativistic wave function [17) . 



2.3 Relativistic MBPT. QED effects. 

Relativistic many-body calculations normally start from the projected Dirac- 
Coulomb-Breit Hamiltonian [T5|, using the Coulomb gauge. 



N 

Hnvpa 



iV 

= A+[j2'^D{i) + Vc + VB\A+, (21) 
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Figure 3: Graphical representation of the all-order Coulomb pair function. The "folded" term 
refers to the intermediate model-space contributions, discussed in the text. For a two-electron 
system with no electron core this represents exactly the non-rclativistic correlation effect. 




Figure 4: Effects beyond the NVPA are conventionally referred to as "QED effects". The first 
diagram represents retardation, which together with the effect of virtual electron-positron 
pairs is referred to as non-radiative QED effects. The following diagrams represent radiative 
effects. The second diagram represents the electron self-energy, the third vertex correction, and 
the last three vacuum polarization/photon self-energy. The internal vertical lines represent 
electron propagators, with positive- and negative-energy orbitals. 



where 

"^^^^^^—^ 4 J 

is the Breit interaction, is the Dirac alpha matrix vector for particle i and 
A4. is a projection operator that excludes negative-energy states. This is known 
as the No- Virtual-Pair Approximation (NVPA) . The all-order procedure de- 
scribed above can also be applied in this scheme. Effects beyond NVPA are 
conventionally referred to as QED effects (see Fig. |4|. 

3 Energy-dependent perturbation theory 

Next, we turn to time-dependent perturbation expansion in the form we have 
developed [H [H [101 [HI [El US] . The standard time-evolution operator de- 
scribes the time evolution of the Schrodinger wave function in the interaction 
picture (IP) 

Mt) = Uiit,to)Mto)- (23) 
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Particles 1 





Figure 5; Comparison between the standard evolution operator, the Green's function and 
the covariant evolution operator for single-photon exchange in the equal-time approximation. 
(Prom [3 Fig. 6.1]). 



The IP is defined by the transfornition 

^i{t)^e'"°*/^^{t), (24) 

where represents the standard Schrodinger representation. The evolution 
operator satisfies the equation 

ih^Ui{t,to)^ViUi{t,to), (25) 

where Vi is the perturbation in IP. 

The standard evolution operator is relativistically non- covariant, since the 
time can flow only in the positive direction. It is for a two-particle system with 
a single photon exchange represented by the first diagram in Fig. |5] 

The single-particle Green's function can for a one-dimensional model space 
be defined 

^ ' (Oh|Oh) ^ ' 

where T is the Wick time-ordering operator and |0h) is the vacuum state, and 
ipii{x), '4'li{x) are the electron-field operators in the Heisenberg representation 
(HP), defined by 

■q,^(t) = e'^*/'"'*(i). (27) 



The numerator of Eq. ( 26 ) is normally singular, but the singularities are in 
this case eliminated by the denominator, so the ratio is regular. 

The Green's function is relativistically covariant . It is for a single-photon 
exchange between two electrons represented by the second diagram in Fig. [Sj 
The free lines are represented by electron propagators, where time can flow in 
both directions. 

We can define a covariant evolution operator (CEO) by generalizing the 
standard evolution operator, as illustrated for single-photon exchange in the last 
diagram in Fig. [5] Here, electron-operator lines are inserted at the free ends 
of the Green's-function diagram, transforming this into an operator diagram. 
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The single-particle CEO (in IP) can be generally defined in analogy with the 
Green's function Eq. (26 1 as 



This can be singular, like the numerator of the Green's- functions definition (26). 
Since we want the procedure here to be valid also for a general multi-dimensional 
model space, we shall use a different procedure for the regularization. 

We define a Green's operator (GO), Q{t,tQ), by the relation (in the following 
omitting the subscript 'Gov') 

U{t,to)P^git,to)-PU{0,to)P. (29) 

Here, the heavy dot implies that the denominators to the left of the dot are 
related to the model space at the position of the dot. This is analogous to the 
procedure we used in defining the counterterm in the time-independent case 
(15). It can be shown that the GO is regular at all times. 

We assume that we have an adiabatic damping, so that at the limit t — )• — oo 
the wave fundtion becomes an eigenfunction of the model Hamiltonian ([T]), 

lim ^'"(i) = (30) 

a state we refer to as the parent state. According to the Gell-Mann-Low theorem 
(GML) [inilin], the wave function at t = can be expressed as the limit 

when the damping vanishes. Here, the numerator and the denominator are both 
normally singular, and only the ratio is regular. 

We assume, in analogy with the non-relativistic case, that the CEO describes 
the evolution of the relativistic wave function 

^{t) = Uit,to)^{to), (32) 

so according to the Gell-Mann-Low theorem 

*(t) =iV„C/(i,-oo)$„, (33) 

where Na is the normalization constant 

1 

(vI/E?|f/(0,-oo)|$")- 



According to the definition ( p9| we have 

C/(0,-oo)P = g(0,-oo)-P[/(0,-oo)P, (35) 
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so the wave function at i = becomes 

^-"(0) = 7V„^;(0,-oo).PC/(0,-oo)$a. (36) 

But NaPU{0, — oo)$Q is according to the GML theorem equal to the model 
function ([5|, which leads to the important relation 

*"(0) =^(0,-oo)*^. (37) 

Comparing with the relation Q, we then find that ^(0, — oo) is the analogue of 
the MBPT wave operator 

n = 5(0, -oo). (38) 



We shall now look at the perturbation expansion of the Green's operator 
( 29 ) . For t = we have C/(°)P = g(^)p = P, which gives from the definition 

g^^'>P = QU^^'>P, (39) 

g^"^^ p ^ Qu'^^'i p - g^-^^ ■ pu^^^ p, (40) 
g'-^^p = Qu'-^^p - g'-^^ ■ pu'^'^^p - g(^) • pu'-^^p, (41) 

etc. The negative terms are counterterms that eliminate the singularities of the 
evolution operator in close analogy with the treatment of the previous section. 
Using Qt/(2) ^ TqV{Tq +Tp)VP we have 

t/(2)F = t/f + A(c;W)Pc/Wp, 

where Go = ^ + ^qV + ^qV^qV + • • • is the GO without intermediate model- 
space states, and 

A(g(i))pf/(i)F = g'-^H£)P£'U'^^^P£ - g^'^\£')P£'U'^^^P£. 

Similarly, we have 

g(3)p = g(i)Qi7(2)p + A(g(i))PC/(2)p - . p[/(i)p 

We want to express tJ^^^P as g^'^^P = TqVg^'^'^P + X in order to find an 
analogue of the MBPT Bloch equation ( 11 ). This gives 

X = g'-^^pu^^^p - g^'^ • PJ7(1)P + A(g(i)) (pj7(2)p - pc/(i)p . P[/(i)p) - 
5(i)A(^;(i))p[/(i)p 

= A(5f )M/o^i' + A(g«)W^(') + A(A(g(i))M/o^i))Tyo(i) -g(i)A(g(i))M/o^i' 



= A*(5(2))^a) ^^(^a))^a) 

where the asterisk indicates that the difference is taken only with respect to the 
last interaction TqV . This leads to 

g(2) = rQ^c;W + ^^T4^« (42) 
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(43) 



s*g 

g = i + rQvg + ^w. 



which is consistent with the general equation fH', Eq. (66)], [7| 6.106] for t = 

(44) 

This can be expressed order by order as 

(45) 



c^Q (m) 



6£ 



This is the Bloch equation for general energy-dependent perturbation. For an 
energy-independent perturbation the differentiation of the last interaction TqV 
leads to a factor of (— Fg), and Eq. (11) is retrieved. If the perturbations can 
be different in different orders, g^"^^ in the sum should be formed by the m last 
interactions and W by the remaining ones, also in complete analogy with the 



MBPT case (111 



4 Combined QED-MBPT 

With the energy-dependent perturbation theory that we have developed, and- 
described in the previous section, we can mix (one-photon) energy-dependent 
QED perturbations with Coulomb interactions to arbitrary order [21] [8] . 

We employ the Coulomb gauge in order to be able to utilize as much as 
possible of the development in standard MBPT. This has among other things 
the advantage that we can represent most of the transverse-photon exchange by 



means of the instantaneous Breit interaction Eq. (22). 

In Fig. [6] we illustrate the incorporation of a transverse-photon exchange, 
i.e., Breit interaction (retarded and unretarded), into an all-order pair func- 
tion, and in Fig. [7| the same when also a Coulomb interaction is crossing the 
Breit interaction. By replacing one of the Coulomb interactions by an instanta- 
neous Breit interaction, we can also take account of the double Breit interaction 
where one interaction is instantaneous and one is retarded. Furthermore, mul- 
tiple Coulomb crossings as well as virtual pairs (negative energy states) can be 
handled. These two sets represent the leading non-radiative QED interactions 
in combination with all-order electron correlation. They have been evaluated for 
He- like ions in the ground state, Is^, by Hedendahl in his PhD thesis [9 . The 
diagrams beyond second order had never been calculated before. Some results 
are shown in Table IH 

To insert the radiative effects into the all-order pair function is more difficult, 
but such work is now in progress at our laboratory, and some preliminary results 
for the self-energy part are inserted in Table|4] As mentioned, it is advantageous 
to use the Coulomb gauge in this type of calculation, but numerical evaluation 
of the self-energy in this gauge has not been done until quite recently. The first 



11 















_____ J 


1 


_ 1 








+ 








1 


1 










_____ H 




_____ H 


1 


_ 1 



Figure 6: Graphical representation of tlie all-order pair function combined with a transverse 
(Breit) interaction. The first two diagrams represent one- and two-photon exchange that 
have previously been calculated, while the remaining ones represent higher-order effects not 
calculated before. 




Figure 7: Same as Fig. |6]with a crossing Coulomb interaction. 



Table 4: QED-correlation effects for He-like ions beyond two-photon exchange 
(first four columns from Hedendahl [Hj, in eV). 



z 


Full Breit 


Ret. part 


Cross. Coul 


Virt. Pairs 


Self-energy 


10 


0.0061 


-0.0011 


-0.0006 


0.0002 




14 


0.0082 


-0.0019 


-0.0010 


0.0004 




18 


0.0191 


-0.0027 


-0.0014 


0.0006 


(-0.004) 


24 


0.0127 


-0.0042 


-0.0021 






30 


0.0150 


-0.0057 


-0.0028 


0.0014 




36 










(-0.010) 


42 


0.0187 


-0.0087 


-0.004 


0.0019 




50 


0.0215 


-0.011 


-0.005 


0.0024 




66 


0.025 


-0.015 


-0.006 


0.0030 






Figure 8: Graphical representation of the all-order pair function combined with the electron 
self-energy. 
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Figure 9: Graphical representation of the all-order pair function combined with the Coulomb 
vertex correction. 




Figure 10: Examples of reducible multi-photon diagrams that can be evaluated by iterating 
the procedure with various QED perturbations together with Coulomb interactions. So far, 
only the first two-photon diagram has been calculated. 



calculation of this kind was performed in 2011 for hydrogcnlike ions by Heden- 
dahl and Holniberg [22| . Inserting the self-energy and the vertex correction into 
the Coulomb pair function, leads to the sequences illustrated in Figs [8] and [9] 

It is possible to insert several (single-photon) QED perturbations in the 
expansion, provided they are reducible, implying that a separated horizontal 
time line can be inserted between the individual interactions. A few examples 
of such combinations are given in Fig. |10[ Again, effects beyond the first two- 
photon diagram have never been accounted for in the past. 

As mentioned, in using the Coulomb gauge, one or several of the Coulomb 
interactions can be replaced by instantaneous Breit interactions Eq. (22 1, and 




Figure 11: Irreducible retarded multi-photon effects can for the time being not be included 
for computational reasons. Most of the effect can be included, though, by replacing one of 
the reatarded Breit interactions by the instantaneous one. 
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in this way most of the effect of the multi-photon exchange can be accounted 
for. Therefore, the procedure presented here can take care of all MBPT-QED 
effects, except irreducible ones, where more than one interaction is retarded, 
as illustrated in Fig. [Tl] These parts, however, which are inaccessible for the 
moment for computational reasons, are expected to be very small, representing, 
maybe, one or a few percent of the total two-photon effect for medium-heavy 
ions and much smaller than some of the higher-order effects in Figs [6] to [9] that 
the procedure does include. Furthermore, the pure two-photon part — without 
electron correlation — can here be evaluated by means of standard two-photon 
programs and added separately [231 IS] • 

We have here considered only two-electron systems, but the procedure can 
be extended to many-electron systems in exactly the same way as in standard 
many-body theory. In addition, it can be combined with the effective Coupled- 
Cluster Approach, as recently demonstrated (CCA) [24]. 

5 Conclusions 

An accurate procedure is now being developed for calculating energy levels and 
separations on few-electron highly charged ions. The procedure utilizes the 
newly developed energy-dependent perturbation theory and can thereby com- 
bine different effects of quantum-electrodynamics (QED) with electron correla- 
tion to all orders in a rigorous manner. The procedure will include important 
effects that have never been calculated before and is therefore expected to yield 
considerably higher accuracy than has previously been attainable. The proce- 
dure is currently being implemented and tested, to begin with on the ground- 
state of He-like ions. The purpose is to test the theory of QED for highly 
charged, few-electron ions on a higher level than has previously been possible. 
To what extent this might explain the discrepancies found by Chantler (Fig. [I]) 
remains to be seen. 
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